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ABSTRACT 

We present some astro-physical consequences of the metric for a point mass in an 
expanding universe derived in Nandra, Lasenby & Hobson, and of the associated in- 
variant expression for the force required to keep a test particle at rest relative to the 
central mass. We focus on the effect of an expanding universe on massive objects on the 
scale of galaxies and clusters. Using Newtonian and general-relativistic approaches, we 
identify two important time-dependent physical radii for such objects when the cosmo- 
logical expansion is accelerating; these radii are found to be insensitive to relativistic 
effects. The first radius, rp, is that at which the total radial force on a test particle 
is zero, which is also the radius of the largest possible circular orbit about the central 
mass m and where the gas pressure and its gradient vanish. The second radius, r<j, is 
that of the largest possible stable circular orbit, which we interpret as the theoretical 
maximum size for an object of mass m. The radius r$ is typically smaller than tf by 
a factor ~ 1.6. In contrast, for a decelerating cosmological expansion, no such finite 
radii exist. Assuming a cosmological expansion consistent with a ACDM concordance 
model, at the present epoch we find that these radii put a sensible constraint on the 
typical sizes of both galaxies and clusters at low redshift. For galaxies, we also find 
that these radii agree closely with zeroes in the radial velocity field in the neighbour- 
hood of nearby galaxies, as inferred by Peirani & Pacheco from recent observations 
of stellar velocities. We then consider the future effect on massive objects of an accel- 
erating cosmological expansion driven by phantom energy, for which the universe is 
predicted to end in a 'Big Rip' at a finite time in the future at which the scale factor 
and the Hubble parameter become singular. In particular, we present a novel way of 
calculating the time prior to the Big Rip that an object of a given mass and size will 
become gravitationally unbound. 

Key words: gravitation - cosmology: theory - galaxies: evolution - galaxies: kine- 
matics and dynamics - galaxies: general - galaxies: clusters: general 



1 INTRODUCTION 



In Nandra, Lasenby & Hobson (2011 1 (hereafter NLH1), we 
considered the effect of a massive object on an expanding 
universe by deriving the metrics that describe a point mass 
in spatially-flat and spatially-curved cosmologies. In partic- 
ular, we derived a general invariant expression for the force 
required to keep a test particle at rest relative to the point 
mass. Conversely, one would expect cosmological expansion 
to have a significant effect on the dynamics of massive ob- 
jects. In this paper, we investigate this effect for objects on 
the scale of galaxies and clusters. 

Considering the radial motion of a test particle in a 
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spatially-flat expanding universe, we showed in NLH1 that, 
in the Newtonian limit, the radial force F per unit mass at 
a distance r from a point mass m is given by 



F = -y 2 -q{i)H\t)r, (1) 

where H(t) = R'(t)/R(t) is the Hubble parameter, q(t) = 
—H'(t)/H 2 (t) — 1 is the deceleration parameter, R(t) is the 
universal scale factor and primes denote derivatives with re- 
spect to the cosmic time irjThus, the force consists of the 
usual time-independent 1/r component due to the central 
mass and a time-dependent cosmological component pro- 
portional to r that is directed outwards (inwards) when 
the cosmological expansion is accelerating (decelerating); see 



also Davis, Lineweaver & Webb ( |2003[ ). For the currently 



(MPH) 



1 We will work in natural units throughout, so that c = G = 1. 
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favoured spatially-flat ACDM cosmology, the cosmological 
force reverses direction at redshift z « 0.67. 

We will consider this general case (and its general- 
relativistic extension) in due course, but the main physical 
effects of the present-day accelerating cosmological expan- 
sion on the dynamics of massive objects can be illustrated 
by considering the special case of a de Sitter universe, which 
contains no cosmological matter fluid (or radiation), but 
only dark energy in the form of a non-zero cosmological con- 
stant A. In this case, the Hubble parameter and, hence, the 
deceleration parameter become time-independent and are 
given by H = \/A/3 and q = — 1. Thus, the force (JlJ) also 
becomes time-independent and reads 



F 



171 i . 

72 + 3 Ar - 



(2) 



The time-independence of this expression is related to 
the well-known (but curious) fact that the presence of a pure 
non-zero cosmological constant can be treated mathemati- 
cally in a completely static manner. As we will discuss in 
Section [4j even in the full general-relativistic case, the met- 
ric for a point mass embedded in an expanding de Sitter 
cosmological background is connected by a simple coordi- 
nate transformation to the standard Schwarzschild-de Sit- 
ter metric, which is static. Thus, in this special case, this 
correspondence unburdens us from the need to include di- 
rectly the expanding cosmological background and allows 
us to consider simply the static case of a point mass in 
the presence of a non-zero cosmological constant. Indeed, 
|2]) is easily obtained by considering directly the Newto- 
nian limit of the Einstein field equations with A / in 
this case (McCauley 19971. Moreover, this correspondence 



enables us to generalise straightforwardly to spatially finite 
(i.e. non-pointlike) spherically-symmetric massive objects. 
In the Newtonian limit, for example, |2| is replaced simply 
by 



F 



„2 ^ S-*"' 



(3) 



where M(r) is the total mass of the object contained within 
the radius r. If the object has the radial density profile p{r ), 
then M(r) = JJ" 4irf 2 p(f) dr. 

Although the de Sitter background is not an accurate 
representation of our universe, the standard cosmological 
model is dominated by dark-energy in a form consistent with 
a simple cosmological constant, which is driving the current 
observed acceleration of the universal expansion ( Straumann| 
2002l|Hubble|1929||Riess et aJ.|2004||Kratochvil et al.|2004| 
Garnavich et al.|1998[ ) . It is therefore of interest to consider 
the astrophysical (as opposed to cosmological) consequences 
of a non-zero A, which has previously been considered only 



in the context of growth of structure analyses ( |Lahav et al. 
[T99lj ). 

Even in the simple Newtonian case, embodied by the 
force expression we see immediately that there is an 
obvious, but profound, difference between the cases A = 
and A 7^ 0. In the former, the force on a constituent par- 
ticle of a galaxy or cluster (say) is attractive for all values 
of r and tends gradually to zero as r — > 00 (for any sensi- 
ble radial density profile). In the latter case, however, the 
force on a constituent particle (or equivalently its radial ac- 
celeration) vanishes at the finite radius vf which satisfies 
tf = [3M(7\f)/A] 1,/3 , beyond which the net force becomes 




Figure 1. Radial-velocity versus distance of stars in the M81 
group, which has a mass of 9.2 X 10 11 Mq, and the best-fit line 
^ (reproduced from |Peirani fc Pachcco 2008). 



repulsive. This suggests that a non-zero A should set a max- 
imum size, dependent on mass, for galaxies and clusters. 

Some observational evidence for this suggestion has 
been reported by Peirani & Pachcco (2005, 2008), who anal- 
ysed the radial velocities of stars in the neighbourhood of 
a number of nearby galaxies, including M83, M81, IC 342 
and NGC 253, by collating a number of unrelated obser- 
vational studies (Karachcntsev 2005]). Note that all these 
galaxies lie at very low redshift (z ~ 10 -4 -10 -3 ). The data 
for M81 are shown in Fig. [I] To obtain the best-fit line for 
each galaxy, it was modelled as a system consisting of low- 
mass satellites sitting within expanding shells around a dom- 
inant core of mass m. Such an approach had already been 
proposed by |Lynden-Bell| fT98"T} and |Sandage| (|1986[) and is 
based on the Lemaitre-Tolman (LT) model ( |Lemaitre|1933 



Tolman|1934 1. One essentially begins with the simple New- 
tonian expression Q for the radial acceleration from which 
an expression for the radial velocity can be calculated, con- 
taining initially unknown constants. Obtaining an analytical 
expression in the case of a non-zero cosmological constant is 
problematic, so Peirani & Pacheco instead found a numer- 
ical solution to the relativistic LT model, valid for z ~ 0: 



Wradial(f) = — 



0.976 / m 



H 2 



3 



1.377ir r, 



(4) 



where n = 0.627 and Ho is the present-day Hubble con- 
stant. This relation appears to fit the data extremely well, 
and highlights the existence of a well-defined 'cut-off radius, 
depending on the galaxy mass, at which the radial velocity is 
zero, separating bound and unbound material. We see that 
this radius is ~ 0.8 Mpc for M81. This is around 3 times 
the typical extent of ~ 0.26 Mpc quoted in |Klypin, Zhao| 



mass to M81. 



& Somerville ( 2002 I for the Milky Way, which has a similar 



The numerical analysis of Peirani & Pacheco (2005 



2008 1 is based on radial velocities, and it relies on observa- 



tional data to fit values to the initially unknown constants. 
To avoid this reliance, in this paper we consider instead test 
particles in circular orbits about the central mass, which 
is also simpler and potentially more representative of the 
trajectories of constituent particles of galaxies and clusters. 
From ([3]), in the Newtonian limit, the speed of a particle in 
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a circular orbit of radius r is given simply by 



v(r) = 



M(r) 
r 



: Ar 2 , 



(5) 



from which it is clear that no circular orbit can exist beyond 
the radius rp. A second physically significant radius rs is 
that of the largest stable circular orbit, which lies somewhat 
within rp, as we will discuss in Section [4] Indeed, the latter 
provides a more physically meaningful limit on the maxi- 
mum sizes of galaxies and clusters. 

The existence and stability of circular orbits around a 
central mass in an expanding universe have been considered 
previously by |Sussman fe Hernandez ( 2003 1, |Nowakowski, 



Sanabria fc Garcia| (|2002[), |Balaguera-Anoh'nez, Bohmer fc 



Nowakowski ( |2005[ ) and |Nowakowski fc Balaguera-Anolfnez 



1 2006 ) . Some details of their work are given in Section 



4.4 



In particular, they assumed a de Sitter cosmological back- 
ground and treated the central mass as pointlike. In this 
paper, we will relax both these assumptions, although not 
simultaneously. 

The overall structure of this paper is as follows. In Sec- 
tion]^ we perform a circular velocity analysis based on the 
simple Newtonian expression |5| applied to two assumed ra- 
dial density profiles p(r) for the central object: a decaying 
exponential and an NFW profile (|Navarro, Frenk fc White] 



1996| |Matos, Nunez fc Sussman||2004[ ). We compare our re- 
sults with the radial velocity analysis of |Peirani fc Pacheco| 
( |2005| |2008[ ) and also derive an analytical expression for the 
pressure profile. In Section [3] we extend our analysis to the 
general-relativistic case, for which we obtain an analytical 
circular velocity expression in the linearised case and em- 
ploy a numerical approach when considering the full Einstein 
equations. In the latter case, we again consider the pressure 
profile in the central object. In Section [4j we remove our as- 
sumption of a de Sitter background cosmology and consider 
the currently-favoured ACDM model with f2 m ,o = 0.3 and 
Oa,o = 0.7, albeit at the cost of having to model the central 
object as a point mass. By comparing with the de Sitter 
background model in this case we thus obtain 'correction' 
factors that may be applied to our earlier results for the ra- 
dius 7\f- We also obtain an expression for the radius rs of the 
largest stable circular orbit. In Section [5j we investigate the 
more speculative scenario that the expansion of the universe 
is driven by phantom energy, leading to the prediction that 



the universe will end in a 'Big Rip' ( Caldwell, Kamionkowski 
|fc Weinberg|2009| ). In particular, we use our force expression 
to find the times prior to the Big Rip at which objects of a 
given mass and size would be expected to become unbound. 
Our conclusions are presented in Section [6] 

Finally, we note that in this paper we model galax- 
ies/clusters as being composed of a single 'phenomenolog- 
ical' fluid, with a single overall density profile and a single 
associated (effective) pressure required for stability; this is 
explained in more detail in the next section. This avoids the 
complexity of an explicit non-linear multi-fluid treatment, 
whereby one would separate the fluid into its baryonic and 
dark matter components. The analyses we perform in this 
work are therefore designed to highlight general features, 
and only indicate the existence and approximate locations 
of the radii rp and rs- 



Quantity 


Galaxies 


Clusters 


r (kpc) 


3 


500 


r T (Mpc) 


0.26 


5 


M{r T ) (M ) 


10 12 


10 15 



Table 1. Summary of the values assumed in our modelling of 
galaxies and clusters. Note that the values for galaxies correspond 
to the virial radius/mass of the Milky Way, as given in Table 2 
of |Klypin, Zhao fc Somerville| f2002| ) . 



2 NEWTONIAN DE SITTER ANALYSIS 

We first adopt the simple Newtonian expression (|5| for the 
circular velocity. To plot the circular velocity function we 
must choose a form for the density distribution. The cen- 
tral massive object is modelled as a total 'effective' fluid, 
made up of two components: baryonic gas and dark mat- 
ter. It is assumed that both components follow a com- 
mon density distribution (up to an overall normalisation), 
such that the total density p(r) — p s (r) + pdm(r), where 
Pg(r) = /gPdm(r)/(f - / g ) and f e = p s (r)/p(r) is the uni- 
form gas fraction throughout the object. Specifically, assum- 
ing spherical symmetry, we consider two model radial pro- 
files. The first is a decaying exponential 



p{r) = p e -o . 



(6) 



where ro is a characteristic physical scale of the object and 
po is its central density. The second is a Navarro, Frenk and 



White (NFW) profile (Navarro, Frenk fc White|1996||Matos, 
Nunez fc Sussman|2004 1 



p(r) = 



por 



(7) 



where ro is again a characteristic physical scale, but po is a 
characteristic density in this case, since the central density 
of Q is infinite. In each case, we fix the constant ro di- 
rectly and then determine po by assuming the mass M{tt) 
inside a 'typical' size rx for low-redshift galaxies and clus- 
ters as listed in Table [T] We also assume the value of the 
cosmological constant to be A = 10 -35 s -2 , which is consis- 
tent with observations ( Carmeli & Kuzmcnko 200"Tj |Carroll 



[200l| |Wem~berg|1989| |Krauss|2003| ) 

The resulting Newtonian circular velocity |5} is plotted 
in Fig. [2] As expected, in each case the circular velocity, and 
hence the radial force, falls to zero at a finite radial distance 
from the centre of the object. This is in sharp contrast to 
the situation when A = (also plotted in Fig. [2J , in which 
case there is no definite cut-off radius. Although, for A 7^ 
0, the radial force vanishes at r-F, this marks the position 
of an unstable equilibrium, since any radial displacement 
of a test particle would move it to a region where either 



the gravitational or cosmological force dominates (Faraoni & 
Jacques|2007 |. We note that for the galaxy and the cluster, 
the cut-off radius rp is smaller for the exponential radial 
density profile than for the NFW profile. Table [2] gives the 
cut-off radii tf for the modelled objects. 

One may also use the plots in Fig. [2] to read off the 
maximum circular velocities v m&K , which clearly lie close to 
the centres of the modelled objects. Using a decaying expo- 
nential profile j6j, for a typical galaxy v max ~ 530 km/s, 
and for a typical cluster u ma x ~ 1280 km/s. Since the to- 
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(c) Galaxy with NFW p(r) 



Figure 2. Circular velocity v(r), in km/s, obtained from the 
(dotted line). 
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expression |[5jl assuming A = 10 35 s 2 (solid line) and A = 





Galaxies 


Clusters 


Exponential p(r) 


1.06 


10.5 


NFW p{r) 


1.20 


12.0 



Table 2. The distance rp (in Mpc) at which the total radial force 
and circular velocities vanish for 'galaxies' and 'clusters', using a 
Newtonian de Sitter analysis and the values for parameters given 
in Table ^ Note that the rp values are read from the plots in 
Fig. § 



tal effective fluid is largely dominated by dark matter, an 
NFW profile |7| is likely to be more realistic. This leads to 
Wmax ~ 300 km/s for a typical galaxy and i> max ~ 1120 km/s 
for a typical cluster. Indeed in this case the maximum circu- 
lar velocity for a galaxy is only about 1.3 times higher than 
that found bylKlypin, Zhao & Somervillel 12002b. 



Having laid down two possible density profiles, we may 
find the associated pressure p(r) required to keep the over- 
all fluid stable. This too is made up of two components: 
an ordinary pressure associated with the gas p s (r) and an 
effective pressure associated with the dark matter Pdm( r )- 
The latter arises when the motions of dark matter parti- 
cles have suffered phase-mixing and relaxation, such that 
their velocities are randomised and the object is virialised, 



as explained by Lynden-Bell ( 1967 \ and more recently by 



Binney & Tremaine (2008). The total pressure is then given 
by p(r) =p g (r) +p e dul (r). 

The controlling equation for gas in hydrostatic equilib- 



d$(r) 
dr 



(8) 



where $(r) is the gravitational potential. The correspond- 
ing equation for the dark matter, modelled as composed of 
collisionless particles which have achieved dynamical equi- 
librium, is the Jeans equation, 



d 
dr 



dr 



(9) 



where cr r {r) is the radial velocity dispersion (Tormen, 



Bouchet fc White|1997 1 . Note that we have ignored a possi- 
ble 'softening' correction by assuming an isotropic velocity 
distribution, which follows naturally from our assumption of 
full spherical symmetry. Comparing |9| with qSp , one can see 
that the two equations are in identical form, and therefore 
that the effective pressure of the dark matter is given by 

PdL{r) = Pdm(r)ar(r). (10) 

We assume that the only interaction between the gas 
and the dark matter is gravitational, so $(r) corresponds 
to the total gravitational potential. In the presence of a 
non-zero cosmological constant, for a fixed total mass M(r) 
within a radius r, the gravitational potential is $(r) = 
- / r °° (M(r')/r' 2 )dr' - (l/6)Ar 2 . Thus, adding the two 
equations |8f and Q together leads to the standard New- 
tonian equation for hydrostatic equilibrium, namely 



dp(r) 
dr 



(11) 



For A = 0, we see immediately that dp(r)/dr is always neg- 
ative, tending to zero as r -> oo. For A / 0, however, the 
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behaviour is radically different and using (JsJ we see that the 
pressure profile has a minimum precisely at r = tf- 

For brevity, we will consider only the (possibly more re- 
alistic) NFW profile. With an appropriate choice of bound- 
ary condition, one may in fact obtain an analytical form for 
pir) (albeit rather unwieldy) for general A, which we believe 
has not been reported previously. The boundary condition 
in both the A = and A / cases is that the pressure itself 
vanishes at the point where dp(r)/dr — 0. Thus, for A = 0, 
one has p(r) — > as r — > oo, whereas for A / one requires 
p{tf) = 0. In the latter case, one thus has yet another phys- 
ically important interpretation of the radius r = r>. One 
finds that the pressure profile is given by 



p(r) 



A-K P lrl In 1 + 



r J I2irp r 



1 + 



+ 127rp r dilog ( 1 + — ) + 67rp T 



In 1 + 



To 



2nplrl In 



T() 



2-wplrl 



27rp6Vohi 1 + 



+ 



Sirplrl In 1 + 



+ 271-poTo In [ 1 + — 

T 



1 + - 



ApoTp 

3 r i + 



Jo 



+ c, 



(12) 



l + r ln(t) 
1-i 



dt; it is worth noting that 



\-k 2 as t — > oo. The constant of 

6 



where dilog(l + r) 
dilog(l + r) — > — §[lnr 

integration C is set by the above boundary condition on p(r 
For the A = case, it has the analytical value C = 2-K i r^p^, 
whereas for A 7^ its value was found numerically. 

Rather than the pressure profile p(r) itself, of more in- 
terest, perhaps, is the ratio p(r)/p(r). Assuming a perfect 
gas law and mean molecular weight of m p /2, we can relate 
this numerically to a temperature, in K, via 



p(r) 



p(r) 



2k B 



p(r) 



p{r) 



(13) 



where m p is the proton mass, k B is Boltzmann's constant 
and the ratio p(r)/p(r) may be expressed in Nm/kg (SI 
units), or in natural units which we mostly adopt in this 
paper. This temperature profile is, in fact, that of the gas, 
T(r) = T g (r). This follows immediately from our assump- 
tion that the gas fraction f s is uniform throughout the 
galaxy/cluster, in which case, using (JsJ) and (|9| respectively, 
one has 



P 
P 



Ps 



Pdm 



(14) 



For A = 0, the temperature profile ( 13 1 can be expressed 
in a universal scale-free form as 



T(x) = * GP k ° r ° x mp {(a 



5x 



+ 2 2 (l + 3;) 2 [6dilog(l + 2; 



+ x[ir 2 x 3 + (2tt 



3:c + 1)(1 + :e) ln(l + :c) 

3(ln(l + a;)) 2 -lnT] 
7)a; 2 + (7r 2 -9)a;-l]}, (15) 



where x = r/ro- However, for A 7^ 0, the term in A in 
( 12 1 breaks scale invariance, and we can no longer write the 
temperature profile in this universal way. 



The peak in the temperature profile ( 15 1 occurs at x = 
x m ~ 0.76, and in terms of x m one can show analytically 
that the peak temperature has the value 

x T {I + x T ) [(1 + x m ) ln(l + x m ) - Xm] GM T m p 
max 2x m (l + 3xm) [(1 + it) ln(l + x T ) — x T ] r T k B 

(16) 

where, as above, Mt is the object's mass at radius tt, and 
xt = vt/tq. Since the total gravitational potential energy 
of the object is ~ GA1 t /tt, and the kinetic energy of the 
constituent particles is ~ (Mt /m p )kBT ', it follows that our 
expression for T m 



( 16 1 in fact corresponds to an analyti- 



cal factor times the the virial temperature. For the 'typical' 
values of tt and To for a cluster given in Table [l] xt — 10 
and we obtain 



0.32 



GM T m p 
r T k B 



3.3 x 10 7 K 



(17) 



or T max ~ 2.8 keV. We may compare this calculated value 
against the virial temperatures of similar sized clusters as 
inferred from X-ray data and simulations; indeed our value 
for T max is in close agreement with the simulations of |Ron-| 



carelli et al. ( 2006 1 and lies within the range found in X-ray 
observations by |Vikhlinin et al~| ( |2006[ ). 

Note that for a galaxy with the parameters in Table 1 we 
can calculate that Tinax w 2.4 x 10 6 K. There has been long 
dispute about whether a component of hot halo gas exists 
for galaxies, but recent X-ray observations, summarised in 



Crain et al. ( 2010 1, do support such a component. In fact the 



temperature of the hot gas halo of the Milky Way itself (the 
mass and size parameters of which were used in Table 1) has 
recently been measured l | Henley et al.p OlO ) and lies in the 
range 1.8— 2. 4x 10 6 K, in good agreement with our calculated 
figure. Of course the behaviour of the gas in the inner regions 
of galactic haloes and clusters is more or less independent 
of the effects of A or universe expansion (see Fig. 3, where 
the effects of A only become perceptible to the right of the 
peak in temperature), but we mention these details, and 
corroborating observations, to show that our model does 
work sensibly in these inner regions, and therefore plausibly 
reflects reality in the outer regions as well. 

Turning then to the broader features of the temperature 
profiles, for A = 0, the ratio p(r)/p(r) is plotted in Fig.|3]for 
a galaxy and a cluster (dotted lines), and exhibits an initial 
increase away from the centre, to the peak just described, 
and then a gradual tailing off as r — > 00. This should be com- 
pared with the corresponding temperature profiles shown in 
Fig. [3] for A / (solid lines). We see that the introduc- 
tion of a non-zero A has a similar effect on the temperature 
profile p(r)/p(r) as it had on the circular velocity profile 
v(r), namely it introduces a definite cut-off radius at which 
the temperature (and pressure) vanish. Since both p(r) and 
dp(r)/dr vanish at r = rp, the temperature profile p(r)/p(r) 
also vanishes there. 

Finally, rearranging equation ( 10 1 and employing the 
equivalence of the ratios given by equation (14 1, it is found 



that p(r)/p(r) is also related to the velocity dispersion of 
the dark matter particles through 



a r (r) 



p(r) 



2k B T(r) 



(18) 



where the right-most-side is obtained using (131. Since this 
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(a) Galaxy (b) Cluster 

Figure 3. Temperature profiles p(r)/p(r) for a galaxy and a cluster with NFW density profiles obtained using Newtonian theory assuming 
A = 10 -35 s — 2 (solid line) and A = (dotted line). Note that the solid lines actually rise again beyond their minima, but since these 
results are unphysical we have omitted them in the plots. 



virial relation is independent of the mass of the dark matter 
particles, it is expected to be applicable to all constituents 
of the overall object which behave like collisionless particles 
in the gravitational potential. Using it to calculate the peak 
velocity dispersion of galaxies within a cluster, using the 
peak temperature given by ( 17 1 we obtain <r max m 740 km/s, 



which is of the same o r der of magnitude as the averages 
found by |Fleenor et al" (20061. The overall scale for <r, 2 (r) 
through a whole galaxy /cluster is shown on a second vertical 
axis in Fig. [3] 

Even from this simple Newtonian analysis, in which the 
expanding universe is modelled rather crudely as a de Sitter 
background, comparing Tables [I] and [2] we see that tf ~ 
2rr for clusters and tf ~ 4rr for galaxies, which is very 
encouraging. In addition, for galaxies the values of Tf are 
also in broad agreement with the 'cut-off' radius at which 
the radial velocity field vanishes in the neighbourhood of 



M81, as found observationally by Peirani & Pacheco (2005 
2008[ ) (see Fig. [I). Given the success of our basic Newtonian 
approach, in the next section we investigate the effects of 
general-relativistic corrections on our analysis. 



3 GENERAL-RELATIVISTIC DE SITTER 
ANALYSIS 

We now extend our analysis to include general-relativistic 
effects to obtain a more accurate estimate of the radius tf 
for galaxies and clusters and also to study further the pres- 
sure profile in such objects. We again assume a de Sitter 
background universe. 

Einstein's field equations with the inclusion of a non- 
zero cosmological constant are 



(19) 



where « = 8ir, T M „ is the matter energy-momentum tensor, 
R^ v = R p MP „ and R denote the Ricci tensor and scalar re- 
spectively, and we adopt the metric signature (+, — , — , — ). 
We use two approaches to solve these equations and find an 
expression for the speed v(r) of a test particle in a circu- 
lar orbit of coordinate radius r about a massive object, as 
measured by a stationary observer at that radius. First, we 
linearise the field equations, which necessitates assuming a 
pressureless fluid within the massive object, but yields an 



analytical form for v(r). Second, we consider the full non- 
linear field equations and again include the effect of pres- 
sure; this requires a numerical solution to obtain v(r) and 
the fluid pressure profile p(r) within the massive object. In 
both cases, we restrict the central object to be spherically 
symmetric and again denote M(r) = JT 4ivf 2 p(f)df. 

3.1 Linearised field equations 

The linearised field equations, with their attendant Lorenz 
gauge constraint, are often written in the form 

LTV" = -2(kT"" + Krf"), 

d^ir = o, (20) 

where = ?? M „ + ft M „, rj^ is the Minkowski metric of 
flat spacetime and hp V is a small perturbation ( jHobson^j 
Efstathiou & Lascnby 20061. Also the d'Alembertian op- 
erator is defined by □ = dad", = — ^r/^h and 
h = h%. One must realise, however, that this form repre- 
sents the field equations specifically in Cartesian coordinates 
(x,y,z), whereas we wish to work in spherical polar coor- 
dinates (r,6,cj>). In general one may interpret the linearised 
field equations as describing a rank-2, symmetric 'gravita- 
tional field' hpv in Minkowski spacetime, so to convert them 
(and the Lorentz gauge condition) to curvilinear coordinates 
in the same spacetime one simply makes the replacements 
rjfty — > 7 M „ and <9 M — > throughout (Misner, Thorne & 
Wheeler 1973). Here 7 M „ is the Minkowski metric in curvi- 
linear coordinates and V M is the covariant derivative defined 
in terms of 7 M „, namely 



pi/ 1 fj.p 

L Pfj.* 1, 1 



(21) 



Working in spherical polar coordinates we specifically use 
7 M „ = diag(l, -1, -r 2 , -r 2 sin 2 ff). 

3.1.1 Tetrad method 

We now obtain a solution to the linearised field equa- 



tions (201 for a static, spherically-symmetric, pressureless 



mass distribution M(r) using a tetrad-based approach. This 
method will also prove useful in solving the full non-linear 
field equations and provides consistency with our approach 
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in NLH1. Note that a full explanation behind this method is 
given in NLH1, so we skip straight to the tetrad definitions 
here. 

For a spherically-symmetric system in the linearised 
case, we may assume a form for the tetrads obtained simply 
by perturbing those for the flat Minkowski metric in spher- 
ical polar coordinates by unknown functions a(r), b(r) and 
c(r): 



e ° = l + a(r), 

e 1 1 = l + 6(r), 

e 2 2 = (1 + c(r))/r, 

e 3 3 = (l + c(r))/(rsin0), 



e° = 1/(1 + o(r)), 
e\ = 1/(1 + b(r)), 
r/(l + c(r)), 
r sin 0/(1 + c(r)). 



2 

e 2 
e 3 3 



(22) 



For a test particle with general four-velocity u, the compo- 
nents in the coordinate basis — [i, r, 9, <f>] are related to 
those in the tetrad frame (denoted by hats) by u M = e^u 1 , 
where dots denote derivatives with respect to the particle's 
proper time r. Thus, we see that for a particle (or observer) 
at rest in the tetrad frame, so that u 1 = [1,0,0,0], then 
t = 1 + a{r) and r — 9 — <f) — 0. Therefore our tetrad frame 
defines the local laboratory of an observer at fixed spatial 
coordinates. Moreover, the three spacelike tetrad unit vec- 
tors e-i (i — 1,2,3) lie in the same directions as the spatial 
coordinate basis vectors e M (fj, — 1,2,3). 

The line-element corresponding to the tetrads ( |22[ ) is 
given by 



ds 2 = 



1 


2 

dt 2 - 


1 


2 

dr 2 — 


V 


.1 + o(r)_ 


l + b(r)_ 








.1 + c(r)_ 



dfT. 
(23) 



where dQ 2 = d9 2 + sin 2 9 d(j> 2 . To first-order in the pertur- 
bation functions, we thus have 

h*" = diag[-2o(r), 26(r), 2c(r)/r 2 , 2c(r)/(r 2 sin 2 9)\. (24) 
Substituting this form into the field equations and Lorenz 



gauge constraint ( 20 1 (appropriately modified to spherical 
polar coordinates) and taking the energy-momentum ten- 
sor to be that of a perfect fluid with zero pressure, yields 
the coupled differential equations (momentarily dropping 
the explicit dependencies on r for brevity) 



a 2 - b 2 - 2c 2 + -(ai - h- 2ci) + 16np + 2A = 0, 
r 

2 8 
a 2 - b 2 + 2c 2 + -(ai - h + 2ci) + -J (6 - c) - 2A = 0, 

2 4 

a 2 + b 2 + -(ai+h) - -a(6 + c) - 2A = 0, 



(25) 



where ai = da(r)/dr, a 2 = d 2 a(r)/dr 2 , and similar for the 
derivatives of b(r) and c(r), and the Lorenz gauge condition 
becomes 



-ai +bi- 2ci + -(b + c) 



0. 



(26) 



a(r), 6(r) and c(r) gives 

a(r) = - \ ' dr + C 2 , 



b(r) = J 8 (r) + Gi, 
, . / 2 (r) C* 3 /i(r) Ci r , , Ar 2 ^ 



where we have defined the integral functions 
h(r) = J[4Tvp(r) - A]r 2 dr 

h{r) = / [2rh(r) + 2rd + 4vrp(r)r 4 + 3Ar 4 ] dr 



(27) 



Hr)=l I ^±^dr, (28) 



and the d (i = 1, 2, 3, 4) are constants of integration. 



3.1.2 Circular motion in the equatorial plane 



We now derive an expression in terms of the tetrads ( 22 1 
for the speed v(r) of a test particle in a circular orbit of 
coordinate radius r about the central mass, as measured by 
an observer at rest at that radius (i.e. the observer defined 
by our tetrad frame). 

In particular, we consider a test particle moving along 
a circular path (r = 0) in the equatorial plane (9 = n/2). 
For such a particle, the components of the four-velocity in 
the tetrad frame may be written 



ii = [cosh tp(r), 0, 0, sinh ip(r)], 



(29) 



where V , ( T ) is the particle's rapidity in this frame, which 
may, in general, be a function of the particle's proper time 
t. The particle's 3-speed, or circular velocity, v in this frame 
is simply v = u 3 /u° = tanh-!/>. 

One may obtain an expression for v in terms of the 
tetrads ( ]22| most simply by using the geodesic equations in 
the tetrad frame directly, which read (see NLH1) 



U + LU j k U U — 0, 

where the spin-connection ui l jh is given by 

t^ijk 2 i.Cijk ~\~ Cjki Cfeij); 
C ij — 6j \Ofj,e v CVC fj, J . 



(30) 



(31) 



Using the tetrads (22 1, the non-zero elements of the spin- 
connection are w 1 oo = OJ°io and cj 1 22 = '^ 1 33 = — c<j 2 12 = 
— oj 3 i3, where 

1 1 + b(r) da(r) 

1 + a(r) dr 
1 l + b(r) ( dc(r) 

r(l + c(r)) \ dr 



(32) 



Solving these equations for the three unknown functions 



Substituting these expressions into (301 and using (291, 
one finds that the geodesic equation for i = 2 is satis- 
fied identically and the equations for i = 0, 3 both reduce 
to tp = 0, which indicates that, as expected, the particle 
has a constant rapidity, and hence constant 3-speed, in the 
tetrad frame. Finally, the geodesic equation for i = 1 gives 
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Exponential p(r) 


NFW p(r) 


Ci 




-47i-p rg(lnr + 1) 


c. 


-47rporg 


\\rl - 47rp r-2 


c 3 








c 4 


47rporg 


^jAr 2 + 4vrporg 



Table 3. Values of the constants in the functions a(r), b(r) and 
c(r) in (|27[l for each considered density profile. 



tanh 2 ?/> = — 0/00/0/33, from which we immediately obtain 



v(r) 



"(1 + c(r)) 



da(r) 
dr 



(l + a(r))( 



dc(r) 
dr 



c(r) 



(33) 



5. J. 3 Results for example radial density profiles 

Since the circular velocity is expected to behave sensibly 
everywhere, the functions a(r), b(r) and c(r) should not di- 
verge at the origin. By considering how their Taylor expan- 
sions about r — may satisfy this requirement, it is not diffi- 
cult to solve for the constants C; (i = 1, 2, 3, 4) for a specific 
radial density profile p(r). Expressions for the constants for 
the exponential profile Q and the NFW profile Q are given 
in Table [3] The resulting circular velocity v(r) profiles are 
practically identical to those plotted in Fig. [2] and so are not 
presented separately. This supports our previous results and 
indicates that linear general-relativistic corrections to the 
(circular) orbits of test particles are unimportant for mas- 
sive objects with densities typical of galaxies and clusters, 
as might be expected. Nonetheless, the analytical results we 
have derived may prove useful when general-relativistic ef- 
fects for astronomical extended density distributions become 
measurable. 



3.2 Full non-linear field equations 

We now consider the full non-linear Einstein equations and 
model the massive object as a perfect fluid with some (non- 
zero) radial pressure profile p(r), as we did for the New- 
tonian case. To solve the field equations, we again use the 
tetrad-based method described in NLHf and employed in 
the previous section. Fortunately, our task is somewhat sim- 
pler for the full field equations than for the linearised equa- 
tions, since we already considered general, time-dependent, 
spherically-symmetric systems in this case in NLHf (see also 



Lasenby et al. ( f998 1), which we now summarise. 

We have considered a general metric element in terms 
of three unknown functions g\{r,i), g2{r,i) and fi(r,t), 



ds 2 = 



9i 



92 



dt 2 + 



2g 2 



dr dt - 



9i 



-r 2 dQ 2 . 



(34) 



We have then shown that, assuming the matter to be a 
perfect fluid with density p(r,t) and p(r,t), Einstein's field 
equations and the Bianchi identities yield the following re- 



lationships between the unknown quantities: 

G 



Lrfl — 


-Gh => h 


= exp {~/ 




p , M 




L r gi = 


F ^+^2~ 


^Ar - 47trp, 


L t g2 = 


M 

G 9^^ + 


|Ar — 47rrp, 


L t M = 


—4ivg2r p, 




L t p = 




) (fi+p), 


L r M = 


47rgi r 2 p, 




L r p = 


-G(p + p). 





- dr 



Here we have defined two linear differential operators 

L t = fidt + Q2d r , 
L r = gid r , 



(35) 



(36) 



and the function F(r,t), the radial acceleration, G(r,t), and 
the intrinsic mass (or energy) interior to r, M(r,t), by 



Ltgi 
L r g2 
M 



G92, 
Fgi, 

1/2 

or (92 



9i+l 



\Kr 2 



(37) 



where A is the cosmological constant. 

In this work, rather than considering a point mass as in 
NLHf, we are considering a static matter distribution (but 
with a non-zero cosmological constant). Thus all functions 
depend only on the radial coordinate r, and the operator Lt 
reduces simply to L t = gz(r)d r . For a given radial density 
profile p(r), one has M(r) = J Q r 47rf 2 p(f) dr , as previously, 
and our goal is to determine the unknown functions fi(r), 
9i(r), g 2 (r), F(r), G(r) and p(r). 

Substituting the above expression for M(r) into the 
LtM equation in (351, one finds g2(r)p(r) — —gi{r)p(r). 
Since the density and pressure must be positive in a realis- 
tic object, one immediately concludes that giir) = 0. It is 
worth noting at this point that only the 'diagonal' tetrads 
ei (no sum on i) defined in NLHf and their inverses are non- 
zero. Thus, the general form of the line-element has reduced 
to the simple form 



1 



A 2 W 



9, 2 {r) 



-dr 



(38) 



Moreover, for a particle (or observer) at rest in the tetrad 
frame, so that its four-velocity components are [u 1 ] — 
[1,0,0,0], we thus have t = fi(r) and r = = <fi = 0. 
Therefore, once again, our tetrad frame defines the local 
laboratory of an observer at fixed spatial coordinates and 
the three spacelike tetrad unit vectors e.t (i = 1,2,3) lie in 
the same directions as the spatial coordinate basis vectors 
e M (jj, = 1,2,3). 

Since 02 (r) = 0, the M and L r g2 equations in (371, and 



the Ltg2 equation in (351 yield, respectively. 



9i{r) 
F(r) 

G(r) 



9i(r) 



— - ^Ar + 4Trrp{r) 



(39) 
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Moreover, once gi(r) and G(r) have been determined, one 
may immediately obtain fi(r) from the L r f\ equation in 
(351. Thus the solution is fully specified by the three non- 
zero functions gi(r), G(r) and p(r). 

We already have an explicit expression for gi (r) in terms 
of known quantities, and an explicit expression for G(r) that 
depends on the unknown function p(r). The pressure profile 
must satisfy the L r p equation in (351, which reads 



dp(r) 
dr 



G(r) 
9i(r) 



[p(r) +p(r)]. 



(40) 



Combining this equation with the G(r) equation in ( 39 1, one 
obtains the Oppenheimer-Volkov equation with a cosmolog- 



ical constant ( Winter||2000 1, 

dp(r) _ (p(r) + p(r))(M(r) + 47rr 3 p(r) - §Ar 3 ) 



dr 



■ (r - 2M(r) - |Ar 3 ) 



(41) 



which directly relates the pressure and density profiles of a 
static, spherically-symmetric perfect fluid distribution. The 
Oppenheimer-Volkov equation may be considered as the rel- 
ativistic generalisation of the Newtonian equation of hydro- 



static equilibrium (111, to which it reduces in the appropri- 



ate limits. For a given p(r) and suitable boundary condition 
on the pressure, we may solve (41 1 to obtain p(r), and hence 
G(r). 

It is worth noting that in the special case of a point 
mass, such that M(r) = m = constant, p(r) = S(r) and 
p(r) = 0, one quickly finds 



G{r) 
Mr) 



3 - lAr 

r z 3 



1 - 212 - lAr 2 

r 3 



2m _ 1 Ar 2 



(42) 



The corresponding line-element ( 38 1 in this case reads 



,/s- = f 1 - ^ - |Ar 2 ) <ir 



2m 



f Ar 2 dr 2 - r 2 dQ. 2 , (43) 



which is the Schwarzschild-de Sitter metric, as expected. 



3.2.1 Boundary condition on the pressure profile 

The boundary condition on p(r) required to solve the 
Oppenheimer-Volkov equation (41 1 may be calculated to 



Taylor expansion p(r) — ~^2^' =0 Pnr n about r = 0. For the 
NFW density profile |7f , however, more care is required in 
determining the appropriate expansion, since p(r) oc 1/r as 
r — > 0. Therefore, in this limit, M(r) oc r 2 and from (44 1 we 
have p(r) oc lnr. This logarithmic divergence at the origin 
must now be taken into account by instead using a power 
series of the form 



P( r ) = Yl 



Pn + Pn In 



To 



(45) 



For each density profile p(r), the coefficients in the expan- 
sion for p(r) can be found by substituting the corresponding 



series into the Oppenheimer-Volkov equation (41 1. These co- 



efficients are then used to set the boundar y co ndition on the 
pressure at r = 10 _3 ro, from which point (41 1 is numerically 
integrated outwards to obtain p(r). As in the Newtonian 
case, we impose the boundary condition that the pressure 
itself vanishes where dp(r)/dr — 0. 



3.2.2 Circular motion in the equatorial plane 

We again consider a test particle moving along a circular 
path (r = 0) in the equatorial plane (8 — it/ '2), for which the 
components of the four-velocity in the tetrad frame may be 



written as in ( 29 1. Once again, one may obtain an expression 
for the 3-speed v — tanhi/> of the test particle in the tetrad 



frame by using the geodesic equations (301. In this case, 



the non-zero elements of the spin-connection (31 1 are again 



Woo = ^°io and cj 1 2 2 = ^*33 = 



-U) 12 = —LO 13 



where 



w 00 = 



Ul 22 = 



ffi(r) d/i(r) 
9i(r) 



dr 



= G(r) 



(46) 



where, for w oo> in the second equality we have used the 
L r fi equation in ( |35[ |. 

We again find that the resulting geodesic equation for 
i — 2 is satisfied identically, that the equations for i = 0, 3 
both reduce to ip — 0, and that the equation for i = 1 gives 
tanh 2 t/j — — oj 1 oo/^ 1 33. Thus the speed v(r) of a particle 
in a circular orbit at coordinate radius r, as measured by a 
stationary observer at that radius is 



v(r) 



rG(r) 
9i(r) 



(47) 



3.2.3 Results for example density profiles 

The resulting circular velocity profiles v(r) and pressure pro- 
files p(r) are again indistinguishable from those plotted in 
Figs [2] and |3] using the Newtonian theory, and are therefore 
not plotted separately. Thus, even full general-relativistic 
corrections to are unimportant for massive objects with den- 
sities typical of galaxies and clusters ( Hwang fc Noh|2006 1 . 

Indeed, quite remarkably, one can show that the loca- 
tion of the radius r = tf at which the circular velocity v(r) 
and the pressure p(r) vanish is identical in the Newtonian 
and full general-relativistic cases, even for large densities 
and pressures at which general-relativistic effects become 
as r — > 0, and hence M(r) oc r 3 . From ( 44 1 , we thus see important. From (47l, we see that the radius r = tf is also 
that p(r) oc r 2 as r — > 0. Thus, one may adopt a simple the point at which the radial acceleration G(r) vanishes, 



a high degree of accuracy by considering a series expansion 
of p(r) about r = 0. For a given density profile p(r), we can 
deduce a suitable form for this expansion by considering the 
Oppenheimer-Volkov equation in its Newtonian limit with 
A set to zero, namely 

dp(r) _ p(r)M(r) 

dr r 2 ' ( ' 

which is the Newtonian equation of hydrostatic equilibrium 
( TTJ) with A = 0. 

For the exponential density profile (J6J L p( r) — > constant 
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changing its direction from inwards for r < tf to outwards 
for r > Tf. From ( 40 1 , this coincides with the radius at 
which dp(r)/dr vanishes, which is also where p(r) vanishes 
according to our boundary condition. Thus, evaluating the 
G(r)-equation in (39 1 at r ■ 



tf 



one finds that 



M(r F ) , , 



(48) 



which is precisely the Newtonian condition defining rp- 

Although rj? is a natural physical radius separating 
bound and unbound material in the galaxy or cluster, it is 
unclear whether r = tf can be interpreted as the maximum 
size of a massive object, since this notion depends more on 
the physical stability of particle orbits, rather than merely 
the point where the radial force on the particle vanishes. 



4 GENERAL- RELATIVISTIC ANALYSIS FOR 
A CENTRAL POINT MASS 

We now address the question of orbit stability and also con- 
sider a more realistic, time-dependent model for the back- 
ground cosmological expansion, albeit at the price of mod- 
elling the central massive object as a point mass. 

In NLH1, we derived the metric for a point mass m 
embedded in an expanding cosmological background, for 
spatially-flat and spatially-curved models. In the spatially- 
flat case, which is a reasonable description of our universe, 
and using non-comoving ('physical') coordinates, the metric 
reads 



2m 2 rr 2 , s 

1 r 2 H 2 (t) 

r 



dt 2 + 2rH(t) 1 



2m 



2m 



dr 2 — r 2 dSl 2 . 



dr dt 
(49) 



Note that this metric is only applicable outside r > 2m. 
Nonetheless it is appropriate to use it for our region of in- 
terest r > m. As shown in NLH1, it is connected by a coor- 
dinate transformation to the corresponding metric derived 
by McVittie (19331. It is worth noting that in the special 



case of a de Sitter background, for which H 2 (t) = A/3, the 
coordinate transformation (Arakida 20101 



dt = dt- 



2 -f §Ar») 



\Kr 2 
3 „„ dr. 



(50) 



converts the metric ( 49 1 into the standard Schwarzschild-de 



Sitter metric 1 43 1 (with t replaced by t), as expected. 



4.1 Radial force on a test particle 

In NLH1, we also obtained an invariant expression for the 
force per unit mass required to keep a test particle at rest 
relative to the central mass m. In the spatially-flat case, this 
reads 



/ = 



^~rH 2 (t) 



2 H 2 (t)) 



1/2 



+ 



vH\t){q{t) + l)J 



1 - 



1- 



3/2 



(51) 



In the special case of a de Sitter background, for which 
H 2 (t) — A/3 and q(t) — —1, this reduces to the well-known 
result 



fas = 



|Ar 



1 _ 2m _ 1 Ar 2 

r 3 



(52) 



As discussed in NLH1, the Newtonian approximation 



to (511 in our region of interest m <Si r <Si 1/H(t), and a 



change of sign, allows us to identify the radial force on a 
unit-mass test particle as simply that given in |l]), namely 



F : 



g(t)H 2 (t)r. 



(53) 



Thus, the force consists of the usual 1/r 2 inwards compo- 
nent due to the central (point) mass m and a cosmological 
component proportional to r that is directed outwards (in- 
wards) when the expansion of the universe is accelerating 
(decelerating) . 

The key difference between the general expression ( 53 1 
and the special case of a de Sitter background is that the cos- 
mological force component in the former is time-dependent. 
Indeed, as shown in section ??, for the standard ACDM con- 
cordance cosmology, the cosmological force term reverses di- 
rection at about z w 0.67, changing from an inwards directed 
force at high redshift (when the universal expansion is de- 
celerating) to an outwards directed force at low redshift (as 
the universal expansion accelerates owing to the dominance 
of the dark energy component). For the standard ACDM 
concordance cosmology, the dominance of the dark energy 
component will increase as the universe continues to expand 
and so the expansion will tend to the de Sitter background 
model considered earlier, for which the cosmological force 
term is time-independent. 

The time-dependence of the cosmological force term in 
the general case means that the radius at which the total ra- 
dial force becomes zero, which marks the cut-off point for the 
existence of circular particle orbits, is also time-dependent. 
For a central point mass, this is given explicitly by 



T F (i) 



m 

'qJtjWW) 



(54) 



provided the universal expansion is accelerating, so that q(t) 
is negative. This also marks the cut-off between gravitation- 
ally bound and unbound material, so it may be interpreted 
as the point separating the 'system' from the rest of the 
universe. Clearly, if the expansion is decelerating (i.e. for 
z > 0.67 in the ACDM concordance model), then the force 
due to the central mass m and the cosmological force are 
both directed inwards and so there is no radius at which the 
total force vanishes, and hence no natural cut-off for the size 
of massive objects. Thus, we see that the time-dependent 
cosmological force term in a general expansion will have a 
profoundly different effect on the formation and structure of 
massive objects as compared with the simple special case of 
a time-indepen dent de Sitter background; in the latter case, 
the expression (541 reduces simply to tf = (3m/A) 1,/3 . 



Nonetheless, at low redshifts, where the dark energy 
component is dominant, we might expect that the values 



of tf obtained using (541 will not differ significantly from 
those obtained assuming a de Sitter background. The dark 
energy density Qa(£) = A/[3H 2 '(t)] has the value Qa, o ~ 0.7 
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Galaxy 


Mass (10 n M©) 


P&P cut-off 


rF(to) 


rs(to) 


M83 


21 


1.10 


1.47 


0.93 


M81 


9.2 


0.75 


1.12 


0.70 


IC 342 


2 


0.51 


0.67 


0.42 


NGC 253 


1.3 


0.40 


0.58 


0.36 



Table 4. Radial distances rp(to) and f"s(*o) i n Mpc from the 
centre of each galaxy, assuming the galaxy to be a point mass, as 
compared to the empirical 'cut-off' radius derived by Peirani & 
Pacheco ('P&P') at which the radial velocity field vanishes; see 
text for details. 



at the current epoch t = to and the deceleration parame- 
ter is currently measured at q(to 
Thus 



-0.55 ( |Krauss||2003 |. 
the 'correction factor' one should apply to our ear- 
lier results in Table [2j derived assuming a de Sitter back- 
ground, is only (— f2A,o/g(£o)) 1//3 ~ 1.08, assuming that the 
correction required in the case of an extended central mass 
is approximately the same as that for a point mass. 

We may compare the present-day values rF(to), calcu- 
lated from (54 1, with the cut-off radii obtained in the nu- 



merical analysis by Peirani & Pacheco ( 2005 2008 1 for spe- 



cific low-redshift galaxies. This comparison is listed in Ta- 
ble [4] from which we see good agreement between the two 
approaches. It should be noted, however, that the rp(to) 
values are typically about a factor ~ 1.4 larger than the 
corresponding radial cut-offs derived by Peirani & Pacheco 
(2005 2008 1. We now address this discrepancy by consid- 



ering the largest stable circular orbits possible about the 
central mass. 



4.2 Largest stable orbits 

We have mentioned previously that the radius r = Tf at 
which the total radial force is zero, and the circular veloc- 
ity v(r) also vanishes, does not necessarily correspond to 
the maximum possible size of the galaxy or cluster. Indeed 
many of the gravitationally-bound particles inside may 
be in unstable circular orbits. A more physically meaning- 
ful 'outer' radius, which one may interpret as the maximum 
size of the object, is that corresponding to the largest stable 
circular orbit, which we denote by rs. 

The radius rs may be determined as the minimum of 
the (time-dependent) effective potential for a test particle 
in orbit about the central mass m ( |Hobson, Efstathiou fc 
Lasenby|2006"} . For the metric {49}, we found in NLH1 that 
the corresponding geodesic equations lead to an equation of 
motion for a test particle given by 



L 2 ( _ 3m\ m . „ 2/ ^ 2r 2 H{t)H'(t) 



rH (t) 

r / r- ( ' 1 - fe-r 2 iP(t) 



rt 



+ 



H'(t)rJl 



1 - '^-r 2 H 2 (t) 
rH'it) 



^1-^(1- Sf-r^f)) 



(55) 



where L is the specific angular momentum of the test par- 
ticle about the central mass m. We recall that dots de- 
note derivatives with respect to the particle's proper time r, 



whereas primes denote derivatives with respect to the cosmic 
time t. 



Expanding ( 55 1 in m and assuming a weak gravitational 
field and low velocities for the test particle (t w 1, r 
one obtains the approximate Newtonian expression 

d 2 r 
dt 2 ' 



q(t)H 2 (t)r + 



0), 



(56) 



This is simply the Newtonian radial force expression ( 53 1 
with the inclusion of a centrifugal term. 

4-2.1 Analytic results 

The effective potential in which the test particle moves is 



given by the first integral of ( 56 1 with respect to r , and so its 
turning points occur at the r-values for which d 2 r/dt 2 = 0, 
namely the solutions of 



q(t)H 2 (t)r 4 +( 



0. 



(57) 



The two real solutions of this quartic both include the factor 
y/256q{t)H 2 (t)L 6 + 27m 4 > 0, thus placing an upper bound 
on L of 



L < 



27m 4 
' 2h&q{t)H 2 {t) 



(58) 



Using the maximum value of L, the radius of the largest 
stable circular orbit is found to be 



rs(t) = 



Aq(t)H 2 (t) 



(59) 



Comparing this expression with that for rj?(t) in (541, 
we see that although both radii are time-dependent, they are 
related by a constant factor. Thus, at any epoch, the radius 
rs(t) of the largest stable circular orbit lies a factor 4 1//3 « 
1.6 inside the radius 7>(t) at which the total radial force on 
a test particle is zero and the circular velocity v(r) vanishes. 
This additional 'correction factor' brings our estimates of 
the maximum possible sizes of individual galaxies into very 
close agreement with the cut-off radii derived by |Peirani fc| 
Pachecoj ( |2005| |2008| , as shown in Table [I] For a galaxy 



of total mass m « 10 Mq, the radius of the largest stable 
circular orbit at the present epoch is rs(to) ~ 0.72 Mpc, and 
for a cluster of total mass m ss 10 Mq we have rs(to) ~ 7.2 
Mpc. These values are consistent with the observed sizes of 
low-redshift galaxies and clusters respectively; only ~ 2.8 
and ~ 1.4 times larger than the typical sizes quoted in Table 
[l] We note that in the special case of a de Sitter background, 
the expression (591 reduces to rs = (3m/(4A)) 1,/3 . 



4-2.2 Numerical results 

As a check on our analytical results, and to include more 
accurately the effect of the radius rs(t) being a function of 
cosmic time, we perform a numerical integration of the New- 
tonian equation of motion {56}. In particular, we consider a 
test particle launched at the present cosmic epoch t — to, at 
an initial radius j-j from the central mass m, such that its 
orbit is initially circular. This corresponds to the boundary 
conditions dr/dt = d 2 r/dt 2 = at t = to, and thus fixes the 
specific angular momentum to be 



L = 



yjmn + q{to)H 2 {t )rf 



(60) 
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Figure 4. Trajectories r(t) (in Mpc versus Gyr) of a test particle launched at the present cosmic epoch t = to, at different initial r adii 
Ti from a central mass m = 1O 15 M0, such that its orbit is initially circular. The ACDM concordance model universal expansion law l |6l| 
is assumed. 



We take m — 1O 15 M0, which corresponds to a typical large 
cluster, and assume a realistic expansion law R(t), corre- 
sponding to the standard spatially-flat concordance model, 



which is given by ( Hobson, Efstathiou fc Lasenby|2006 1 

i 

'(1 



o(t) 



^A,o) 



sinh 



2 H > 



(61) 



where a(t) = R(t)/Ro is the normalised scale factor, Qa,o = 
0.7 and H = 70 km s _1 Mpc" 1 . 

It is found that there is a critical initial radius, n, c ~ 
5.575 Mpc, above which the resulting particle trajectory ex- 
pands with the universe. Physically this means that a test 
particle in an initial orbit of radius > r%, c would ultimately 
be swept out by the universal expansion. The resulting tra- 
jectory r(t) for an initial radius equal to the critical value is 
shown in Fig. |4|a), whereas the trajectories for larger val- 
ues of n are shown in the other panels. One sees that for 
the critical value, the radius of the particle orbit oscillates 
about a fixed value f w 6.3 Mpc that is slightly larger than 
n, which itself lies at the minimum of the oscillation. For 
larger values of ri , the radius of the resulting particle orbit 
grows significantly with cosmic time. Nonetheless, in these 
latter cases, the growth of the particle orbit lags behind 
the universal expansion. Although this 'lag' becomes less 
pronounced for larger values of n, the test particle never 
becomes completely comoving owing to the presence of the 
central mass. 

We note that Price ( 2005 1 has used a very similar 



geodesic equation to (561 to analyse the motion of an elec- 



tron orbiting a central nucleus, with a slight modification to 
account for the dominant electrostatic force in that problem. 
An extension of this work is presented in Appendix ??. 



4.3 Minimum density of objects 

If we consider rs (t) as the maximum possible size of a mas- 
sive object at cosmic time t, and make the simplifying as- 
sumption that objects are spherically-symmetric and have 
constant density, then the condition ( |59[ ) shows that there 
exists a time-dependent minimum density for objects, given 
by 



»(*) = 



3m 



4irr s 3 (t) 



nG 



(62) 



where in the last equality we have reinserted a factor of G 
to convert back to SI units. It should be noted that, since 



(59 1 is valid only for q(t) < (accelerating expansion), then 
this restriction also applies to (621. For q(t) > 0, p m i U (t) 
vanishes. 

We point out that in the special case of a de Sitter 



background, Nowakowski, Sanabria & Garcia (2002) calcu- 



late the minimum density for objects to remain gravitation- 
ally bound to be p m i n = A/(47rG). This is a factor of 4 
smaller than our result for a de Sitter background. It does 
however, correspond to what we would obtain by defining 
Pmin(i) = 3m/ (Airrp ), with rp = (3m/A) 1//3 . Since rp cor- 
responds only to the theoretical maximum size of an ob- 
ject, and not the true maximum attainable size given by 
rs = (3m/(4A)) 1//3 , we believe our estimate of the mini- 
mum density to be more realistic. 



Using 1 62 1 one finds that, at the present cosmic epoch, 
Pmin(io) ~ 4 x 10 -26 kg m -3 , which should be compared 
with the present-day critical density /9 cr it,o ~ 9 x 10 -27 kg 
m~ 3 required for the universe to be spatially-flat. Indeed, 
since at any cosmic epoch, p C rit(t) = 3H 2 (t)/(8TrG), one has 
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simply 

Pcrit (t) 



-8q(t), 



(63) 



provided q(t) < 0. Alternatively, since the average matter 
density p m (t) ~ £lm(t)p C rit(t), the minimum fractional den- 
sity contrast <5 m in(i) = [pmin(i) — Pm{t)]/pm(t) of a matter 
fluctuation at cosmic time t is 



l + 



8g(t) 

n m (t) 



(64) 



Inserting the present-day values fi mi o ~ 0.3 and go ~ —0.55, 
one finds <5 m i n (io) ~ 14. 

The above expressions and values are, of course, rather 
approximate. In particular, our rs values have been calcu- 
lated assuming a central mass point rather than an extended 
mass and realistic astrophysical objects do not have uniform 
densities. Nonetheless, the above arguments can be taken as 
an indication of the existence and rough order of magnitude 
of such a minimum density. Calculated accurately, this min- 
imum density should correspond to that required for any 
structure formation at all to occur at a given cosmic time t 
during an accelerated expansion phase of the universe. We 
leave this as a topic for future research. 



4.4 Comparison with previous studies 



Balaguera-Anolfnez, Bohmer fc Nowakowski| ( |2005| and 
Nowakowski & Balaguera-Anolmez (2006) have previously 



used the Schwarzschild de Sitter metric ( 43 ) to find expres- 



sions for (what we call) rp and rg in the special case of a 
de Sitter background universe. 

It is worth noting that they were able to use a more tra- 
ditional method to find tf and rs, since the Schwarzschild 
de Sitter metric is static. In such cases, one can combine the 
resulting geodesic equations into an 'energy' equation of the 
form |r 2 + U e g(r) = E, where E is the total energy of the 
particle per unit rest mass and one can read off the effec- 
tive potential U a a(r). For the Schwarzschild de Sitter metric 
( 43 1 , this is given by 



U eB (r) 



-Ar 2 
6 



2r 2 



ml/ 



(65) 



The value of tf is then simply the location of the minimum 
of U e s(r) with L = 0, which immediately leads to an ex- 
pression for rp which agrees precisely with our result (541 



in the special case of a de Sitter background. The value of 
rs is the location of the minimum of U e s{r) with L = L max , 
where L max is maximum value of L for which U e s(r) pos- 
sesses a minimum. In the limits L~^> m and A 1 adopted 
by these authors, the resulting expression for rs agrees with 



our result (59) in the de Sitter case. 



This approach should be contrasted with the derivation 



in NLH1 of the equation of motion (551. The latter is more 



complicated owing to the dependence of the metric ( 49 1 on 



the cosmic time t, which prevents one arriving at an 'energy' 
equation. Nonetheless, we note that in the special case of a 



de Sitter background, the equation of motion (55 1 reduces 
to 



ml, L 



3m 
r 



(66) 



The first integral of this equation then immediately yields 
an 'energy' equation with t he effective potential (|65[ ). 

We further note that Sussman & Hernandezl ( 2003 1 



also obtained the factor rp/rs ~ 1.6 for a de Sitter back- 
ground, but individual values of tf and rs are greater than 
ours by a factor of ~ 10 for galaxies and ~ 3 for clusters. 
Their approach deviates significantly from our work and that 
of Balaguera-Anolmez, Bohmer & Nowakowski (20051 and 



Nowakowski & Balaguera-Anolmez (20061, since they model 



galaxies and clusters as dark matter haloes in equilibrium 
with a static 'A-field', rather than simply as point masses. 
It is possible that, by taking the spatial extent of the cen- 
tral object into account, we may close this small discrepancy 
between our results, but we leave this as a topic for future 
research. We point out that Sussman & Hernandez obtain 
the values of vf and rs by numerically integrating two cou- 
pled differential equations, so no analytical results are given 
in their work. 



5 PHANTOM ENERGY AND THE 'BIG RIP' 

We have thus far considered dark energy only in the form 
of a cosmological constant, for which the equation-of-state 
parameter w = — 1. This is well supported by observations, 
but more generic dark energy models allow for w to differ 
from —1, and even to vary with cosmic epoch. In any case, 
W = — 1 is usually considered as the lower-limit, since a 
dark energy 'fluid' with w < — 1 violates all the relativistic 
energy conditions. Nonetheless, this has not deterred some 
theorists from considering the cosmological consequences of 
a dark energy component with a constant equation-of-state 
parameter w < — 1, which is often termed 'phantom-energy' 
( |Caldwell|2002[ ) . We now consider the effect on massive ob- 
jects of a cosmological expansion driven by phantom energy. 



5.1 The Big Rip 



As shown by Caldwell, Kamionkowski & Weinberg (20091, 



the most profound consequence of phantom energy is that it 
leads to a 'Big Rip' in the far future, at which the universe 
scale factor and Hubble parameter become singular. 

To explain this, for simplicity let us consider a spatially- 
flat model, which is a good approximation to our universe. 
The Friedmann equation governing the time evolution of the 
normalised scale factor a(t), assuming negligible radiation 
energy density, is 



H 2 (t) = H 2 [fi m , a- 3 (t) + (1 - fi n 



,o a 



-3(1+™) 



(*)] , (67) 



where H(t) = a'(t)/a(i), ~ 0.3 is the present-day mat- 
ter density parameter and the current value of the Hubble 
parameter is taken to be Ho w 70 km s _1 Mpc -1 . 

At late times the second term on the right-hand side, 
relating to dark energy, dominates over the first, relating to 
matter. Writing 1 + w as — 11 + w\, the Friedmann equa- 



tion may be written as H (t) = H (1 — f2 m ,o)a 
Solving for the normalised scale factor then gives 



a(f) — r 

_3|1+W| (#0 V 1 " SVo(*0 ~t)+ 3|T^J _ 



-3|l+iu| 



(t). 
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(68) 

One sees immediately that a(t), and hence H(t), diverges at 
u 

2 



a finite time t — t r i P in the future that satisfies 



tri 



to 



m,0 

for example, the time 
before the Big Rip is ~ 22 Gyrs. 



(69) 



If w = — |, for example, the time remaining for the universe 



5.2 Times at which objects become unbound 

Since H(t) becomes infinite at t — t r i P , then so does the cos- 



mological component of the radial force ( 53 1 on a test parti- 
cle in the vicinity of a central mass. Moreover, since q(t) < 
for the accelerating expansion during phantom-energy dom- 
inance, this force is directed ouwards, away from the central 
mass. Thus, by the time of the Big Rip, all massive objects 
will have been ripped apart. 

We now consider the time before the Big Rip that this 
occurs for an object of a mass m and size r. Note that 
we only perform a rough calculation due to our simplify- 
ing assumptions. For example, we assume the mass to be 
concentrated at the centre of the object, as in a standard 
Newtonian framework. We then specifically calculate the 
time at which a particle at the edge of the object becomes 
gravitationally unbound from it; that is when the inwards 
force component of magnitude m/r 2 , due to the central 
mass, is balanced by the outwards component of magnitude 
\q(t)H 2 (t)r\ due to the (accelerating) cosmological expan- 
sion. Note that this is not the same as calculating the time 
at which the point mass itself is ripped apart, which we are 
unable to evaluate due to the singularity in our model at 
r = 2m. Nonetheless, it is the time at which the object ceases 
to exist in its current state and begins to be ripped apart, 
as layers of it gradually peel off as the cosmological force in- 
creasingly dominates over the gravitational force. Since the 
particle is at a distance r >> m, the physical singularity in 
our model at r — 2m does not affect the result. 

At late times, substituting a(t) from ( |68[ ) into the ap- 
propriate form of the Friedmann equation H (t) — Hq(1 — 



,o)a 



-3\l+w\ 



(t), one finds 



H\t) = 



31 1 + w\t 



(70) 



where t = t T i p — t is the time before the Big Rip. It is worth 
noting that this result is independent of both H and fi m ,o- 
Moreover, we see that at late times the deceleration param- 
eter tends to the constant value 



9(t) 



2 



■ 3w\ 



(71) 



Thus the time before the Big Rip at which the total radial 
force < |53[ ) on a particle at the surface of an object of mass 
m and radius r vanishes is given by 



t 



311 + ml 



(72) 



Although this is only an approximate result, based on the 
Newtonian force expression ( 53 1 , it remarkably matches 



Caldwell's approximate expression precisely, despite the lat- 
ter being derived in a rather different manner based on en- 
ergy arguments. This could be merely a coincidence, but to 



Object 


m 




r 


t 


Galaxy cluster 


10 15 M 





5 Mpc 


9.3 Gyr 


Milky Way 


5.8 x 10 11 


M Q 


15 kpc 


69.9 Myr 


Solar System 


M 




150 x 10 9 m 


3.6 months 


Earth 


6 x 10 24 


kg 


6371 x 10 3 m 


25.2 mins 



Table 5. The times i before the Big Rip at which objects with a 
given mass m and size r are expected to be ripped apart, assuming 
a phantom energy equation-of-state parameter w = —3/2. In this 
case, the time until the Big Rip is t r ; p — to 22 Gyr. 



determine the true extent to which our two approaches may 
in fact be parallel requires further thought. We leave this for 
future research. 

For w — —3/2, the times t before the Big Rip that 
specific objects of mass m and size r are expected to rip 
apart are given in Table [5] It is assumed (admittedly, rather 
implausibly) that, as the universe evolves, the objects do not 
accrete any additional mass, that their size remains fixed at 
the present-day value until they are ripped apart, and that 
non-gravitational binding forces are negligible. 



6 CONCLUSIONS 



In Nandra, Lasenby & Hobson ( 2011 I we derived the metric 



for a point mass m embedded in an expanding cosmologi- 
cal background, for both spatially-flat and spatially-curved 
universes. We also derived the corresponding invariant ex- 
pression in each case for the force required to keep a test 
particle at rest relative to the central mass. In this paper, we 
have presented some important astrophysical consequences 
of this work, focussing on the effect of an expanding universe 
on massive objects on the scale of galaxies and clusters. 

For a spatially- flat universe in the Newtonian limit, the 
force on a test particle at a distance r from a point mass 
m consists of the usual time-independent 1/r 2 inwards com- 
ponent due to the central mass and a time-dependent cos- 
mological component proportional to r that is directed out- 
wards (inwards) when the expansion of the universe is ac- 
celerating (decelerating). It is immediately clear, therefore, 
that the effect of an expanding universe on a massive ob- 
ject is profoundly different if the expansion is accelerating 
as opposed to decelerating. In the latter case, the force on a 
constituent particle of a galaxy or cluster is attractive for all 
values of r and tends gradually to zero as r — > oo (for any 
sensible radial density profile), with no definite cut-off. In 
the former case, however, the force on a constituent particle 
(or equivalently its radial acceleration) vanishes at a some 
finite radius 7\f (say), beyond which the net force becomes 
repulsive. This is also the radius of the largest possible cir- 
cular orbit about the central mass. This suggests that an 
accelerating cosmological expansion should set a maximum 
size, dependent on mass, for galaxies and clusters. 

Although our force expression accommodates a realis- 
tic dynamical expanding background cosmology, it comes at 
the price of modelling the central massive object as a point 
mass. In reality, it is more appropriate to model galaxies 
or clusters as extended objects having some density profile. 
We consider this in the special case of a de Sitter back- 
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ground, which simplifies matters considerably. In this case, 
the expansion is always accelerating and the analysis can 
be performed in a entirely static manner with the inclu- 
sion of a non-zero cosmological constant. Although the de 
Sitter background is not an accurate representation of our 
universe, the standard cosmological model is dominated by 
dark-energy in a form consistent with a simple cosmological 
constant, which is driving the current observed accelerated 
phase of expansion. 

We calculate the value of tf in this case for spherically- 
symmetric matter distributions with an exponential and 
NFW radial density profile, respectively, using Newtonian 
theory, linearised general-relativity and the full non-linear 
Einstein field equations. The three methods give indistin- 
guishable results for objects with masses and sizes typical of 
galaxies and clusters, as might be expected. More surprising, 
perhaps, is that one can show analytically that the location 
of r = tf is identical in the Newtonian and full general- 
relativistic cases. For clusters we find ~ 10.5-12.0 Mpc 
and for galaxies we find tf ~ 1.06-1.20 Mpc, which, as 
maximum sizes, are compatible with the sizes of 5 Mpc and 
0.26 Mpc respectively. In addition, for galaxies our values 
for tf agree well with the typical radius at which the radial 
velocity field in the neighbourhood of nearby galaxies van- 
ishes, as inferred by |Peirani fc Pacheco| p005} [2008| from 
recent observations of stellar velocities in such objects. In 
the Newtonian and fully general-relativistic analyses, we in- 
clude the effects of pressure and derive the radial pressure 
profile p(r) in the central object. We find that dp(r)/dr also 
vanishes precisely at r — rp, as does p(r) itself according to 
our boundary condition. In the Newtonian case, we obtain 
analytical results for pressure and temperature profiles in 
the case of an NFW density profile. 

The key difference between a general expression and the 
special case of a de Sitter background is that, in the former, 
the cosmological component of the radial force on a test 
particle is time-dependent. This means that the radius Tf at 
which the total radial force becomes zero, which marks the 
cut-off point for the existence of circular particle orbits, is 
also time-dependent. In particular, for the standard ACDM 
concordance cosmology, there is no finite radius at which the 
total force vanishes, and hence no natural cut-off for the size 
of massive objects, during the decelerating phase prior to 
z « 0.67. Nonetheless, returning to our model of the central 
massive object as a point mass, we find that at low redshifts, 
where the dark energy component is dominant, the values of 
tf obtained do not differ significantly from those obtained 
assuming a de Sitter background. 

For a general cosmological expansion, the radius rs 
of the largest stable circular orbit is also time-dependent, 
but we find that at any cosmic epoch it is simply a factor 
4 1 / 3 ~ 1.6 smaller than tf- It is also perhaps a more ap- 
propriate measure than tf for the maximum possible size 
of an object of mass m. Indeed it is encouraging that we 
find the factors rs/rr at the present time for both galaxies 
and clusters to be close to 1; we obtain ~ 2.8 and ~ 1.4 
respectively. We also obtain an approximate expression, de- 
pendent only on cosmological parameters, for the minimum 
density of an object for it to remain gravitationally bound 
at a given cosmic epoch. 

Finally, we also consider the future effect on massive ob- 
jects of cosmological expansion driven by dark energy with 



an equation-of-state parameter w < —1, termed phantom 
energy. In such a model the universe is predicted to end in 
a 'Big Rip' at a finite time in the future at which the scale 
factor and the Hubble parameter become singular. Using 
our radial force equation, we obtain an expression for the 
time prior to the Big Rip that an object of a given mass and 
size will become gravitationally unbound. Remarkably, this 
is found to agree precisely with the approximate expression 
derived by Caldwell using a rather different approach based 
on energy arguments. 
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APPENDIX A: AN ATOM'S 
'ALL-OR-NOTHING' BEHAVIOUR 



Price ( 2005 1 has analysed the motion of an electron of mass 



m e orbiting a central nucleus using the weak field, low ve- 
locity geodesic equation l ]56| . For this model, he has shown 
that it is appropriate simply to replace the m/r 2 gravita- 



tional force in (561 with C/r — QQ'/(47reom e r ), where 



QQ' is the product of the nuclear and electron charges. Price 
questions whether or not the electron eventually takes part 
in the cosmological expansion; that is, does the electron 
follow a trajectory of constant 'physical' coordinate r (no 
atomic expansion) or of constant comoving coordinate f (full 
cosmological expansion of the atom), or does the electron 
do something 'in-between'? Using the model expansion law 
a(t) — 1 + t' 2 tanh(£), where a(t) — R(t)/R(to), and launch- 
ing the test particle in a momentarily circular orbit with 
r = 1 (arbitrarily) at t = 0, (dr/dt) t =o = {d 2 r / dt 2 ) t=Q = 
and thus C = L 2 , Price finds that there is a critical value of 
the electron's specific orbital angular momentum L « 3.46 
that determines an 'all-or-nothing' behaviour: that is a 
weakly bound atom (L ^ 3.46) will eventually begin to 
stretch in size and participate wholly in cosmological ex- 
pansion ('all' behaviour), whereas a tightly bound atom 
(L > 3.46) will eventually ignore the cosmological expan- 
sion and settle into a circular orbit of fixed 'physical' ra- 
dius ('nothing' behaviour). We have been able to reproduce 
Price's results using his expansion law and initial conditions. 




Figure Al. The cosmological radius f = r/a(t) of an electron 
orbit around a central nucleus for the critical value L = 3.46 of 
the electron's specific orbital angular momentum. 



Fig. |A1| shows the comoving radius of the electron orbit as 
a function of cosmic time t for the critical value L = 3.46, 
which shows an apparent plateau at late times that Price 
interpreted as 'all' behaviour. 

Closer inspection of the 'all' behaviour reveals, however, 
a subtlety that was not addressed previously. The apparent 
plateau in Fig. |A1| is not actually flat at all; rather f con- 
tinues to decrease at late times but at a much slower rate 
than previously. For this value of L, the atom, in fact, ex- 
hibits only a 'something' behaviour. Fig. |A2| shows instead 
the 'physical' radius r of the electron's orbit, together with 
the universal scale factor a(t). It then becomes clear that the 
atom does indeed expand forever, but at a rate that increas- 
ingly lags behind the cosmological expansion for L = 3.46. 
As we decrease L, for even less tightly bound atoms, the 
amount by which the electron's expansion lags behind the 
cosmological expansion decreases. Eventually, for L ^ 0.1, 
the two expand at the same rate and the behaviour changes 
from 'something' to 'all'. 

This paper has been typeset from a TgX/ DTjtX file prepared 
by the author. 
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Figure A2. The 'physical' radius r (bold line) of an electron 
orbit around a central nucleus for different values of the electron's 
specific angular momentum: L = 3.46 (top), L = 2.5 (middle), 
L = 0.1 (bottom). The scale factor a(t) is shown as the dashed 
line. 
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